OPEN 3 ACCESS Freely available online 



PLOS 



ONE 



miRNA Expression in Anaplastic Thyroid Carcinomas 

Aline Hebrant\ Sebastien Floor\ Manuel Saiselet\ Aline Antoniou\ Alice Desbuleux\ 
Berengere Snyers\ Caroline La\ Nicolas de Saint Aubain^, Emmanuelle Leteurtre^'^, Guy Andry^, 
Carine Maenhaut^'^* 

1 Institute of Interdisciplinary Research (IRIBHM), School of IVledlcine, Universite libre de Bruxelles, Campus Erasme, Brussels, Belgium, 2lnstitut Jules Bordet, Bruxelles, 
Belgium, 3 Universite de Lille 2, Faculte de IVledecine, Lille, France, 4CHRU de Lille, Institut de Pathologie, Lille, France, 5 WELBIO, School of Medicine, Universite libre de 
Bruxelles, Campus Erasme, Brussels, Belgium 



Abstract 

Anaplastic thyroid carcinoma (ATC) is tine most lethal form of thyroid neoplasia and represents an end stage of thyroid 
tumor progression. No effective treatment exists so far. In this study, we analyzed the miRNA expression profiles of 1 1 ATC 
by microarrays and their relationship with the mRNA expression profiles of the same 1 1 ATC samples. ATC show distinct 
miRNA expression profiles compared to other less aggressive thyroid tumor types. ATC show 18 commonly deregulated 
miRNA compared to normal thyroid tissue (17 downregulated and 1 upregulated miRNA). First, the analysis of a combined 
approach of the mRNA gene expression and of the bioinformatically predicted mRNA targets of the deregulated miRNA 
suggested a role for these regulations in the epithelial to mesenchymal transition (EIVIT) process in ATC. Second, the direct 
interaction between one of the upregulated mRNA target, the LOX gene which is an EIVIT key player, and a downregulated 
miRNA, the miR-29a, was experimentally validated by a luciferase assay in HEK cell. Third, we confirmed that the ATC tissue 
is composed of about 50% of tumor associated macrophages (TAIVl) and suggested, by taking into account our data and 
published data, their most likely direct or paracrine intercommunication between them and the thyroid tumor cells, 
amplifying the tumor aggressiveness. Finally, we demonstrated by In situ hybridization a specific thyrocyte localization of 3 
of the deregulated miRNA: let-7g, miR-29a and miR-30e and we pointed out the importance of identifying the cell type 
localization before drawing any conclusion on the physiopathological role of a given gene. 
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Introduction 

Thyroid tumors are divided into encapsulated benign tumors 
(autonomous and follicular adenomas) and carcinomas. These 
carcinomas are themselves subdivided into differentiated carcino- 
mas (follicular carcinomas (FTC) or papillary carcinomas (PTC)) 
which may evolve into the very aggressive and dedifferentiated 
anaplastic carcinomas (ATC) [1]. Despite its low frequency (<5% 
of all thyroid carcinomas), ATC contributes to 14—50% of the 
thyroid mortality cancer with a medium survival of 3 to 5 months 
[2,3] . Benefits obtained from chemotherapy and radiation therapy 
remain only marginal and there is, until now, no alternative 
treatment [4,5]. New therapeutic approaches are therefore 
needed. On the other hand, the ATC tissues contain tumor 
associated macrophages (TAM) which may represent up to 50% of 
the nucleated cells and show elongated thin ramified cytoplasmic 
extensions [6]. A relationship between the increased density of 
TAM type M2 and decreased survival was reported in thyroid 
cancer patients [6-9] . A direct correlation between the density of 
the presence of the TAM and the aggressiveness of the ATC has 
been suggested [6-9]. 

To identify the molecular mechanisms involved in tumor 
evolution, we have analyzed the mRNA expression profiles of 59 



thyroid tumors ( 1 1 ATC and 48 PTC) using the Affymetrrx 
microarray technology [10]. We have now integrated into these 
results the miRNA expression profiles of the same 1 1 ATC . 

miRNA are non-coding RNA of approximately 22 nucleotides 
in length which can regulate the expression of more than several 
hundreds of genes following therr binding to the 3 ' untranslated 
region (UTR) of the mRNA target containing corresponding 
antisense sequences [11]. The ~6 nt 5' miRNA 'seed' region is 
the most crucial component for recognizing and binding to target 
sites [12]. miRNA regulate gene expression through translational 
inhibition (accounting for 16%) or mRNA degradation (account- 
ing for 84%) [13-16]. Clearly, as a miRNA acts mostly through 
mRNA degradation [15,16], we expect to see this reflected in the 
mRNA levels of its targets and thus in our microarray results. 

As bioinformatics mRNA target prediction algorithms allow to 
predict for each deregulated miRNA a list of mRNA targets, we 
integrated the bioinformatically predicted mRNA target list to the 
deregulated mRNA list of the same 1 1 ATC obtained previously 
experimentally [10]. In this study, we have used the miRtarget2 
algorithm [17] which is based on several criteria: the complemen- 
tarity of the seed region, the observation that true sites tend to be 
conserved between species, the secondary structure (a target site is 
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not likely to be functional if it is inaccessible to miRNA binding) as 
well as the location of the seed target in 3'UTR (a site in the 
middle of a long UTR is less likely to be functional). 

The analysis of miRNA and mRNA gene expression of the 
same ATC samples and bioinformatic target prediction programs 
has highlighted a major deregulated transcription program that 
dictates ATC cells characteristics: the epithelial to mesenchymal 
transition (EMT). Moreover, a luciferase assay validated the direct 
interaction between the 3'UTR of LOX, a key gene in the EMT 
process, upregulated in the 1 1 ATC samples, and miR-29a, 
downregulated in the same 1 1 ATC. 

Because the ATC tissue is composed of about 50% of TAM 
which are suggested to play a major role in the ATC 
aggressiveness, we have confirmed by in situ hybridization the 
thyrocyte derived tumor cell localization of 3 of the deregulated 
miRNA: let-7g, miR-29a and miR-30e. 

Materials and Methods 

Tissue samples 

14 ATC tissues and 19 other thyroid tumors (5 PTC, 5 FTC, 4 
FTA and 5 AA) were provided by different hospitals: Regional 
Reference Cancer Center of Lille (Lille, France), Pitie-Salpetriere 
(Paris, France), Bordet Institute (Brussels, Belgium), Cancer Center 
of the Cliniques Universitaires Saint-Luc (UCL, Brussels, Belgium) 
and Catholic University of Leuven (KUL, Leuven, Belgium). AU 
these tumors were compared to a common miRNA reference pool 
of 19 normal thyroid tissues which are adjacent to benign thyroid 
tumors. 1 1 of these 1 4 ATC samples were used for microarray 
studies and 3 of them wci'c used for independent validation. 
Tissues were immediately dissected, placed on ice, snap-frozen in 
liquid nitrogen and stored at — 80°C until processing. 4 of the 
studied ATC were also used for in situ hybridization and 
immunohistochemistry. Protocols have been approved by the 
ethics committees of ALL institutions. Written informed consent 
was obtained from all participants involved in the study. 

RNA purification 

Total RNA was extracted from thyroid tissues using a Trizol 

reagent kit (Invitrogen), followed by purification on miRNeasy 
columns (Qiagen). The RNA concentration was spectrophotomet- 
rically quantified, and its integrity was verified by visualization 
using an automated electrophoresis system (Experion, Bio-Rad). 

Mutation screening 

In order to determine the mutational status for TP53, BRAF, 
H-RAS, N-RAS, K-RAS, PIK3CA, P-catenin in the 11 ATC 
samples, the mRNA sequences containing the most frequent 
mutations were reverse-transcribed and amplified by PCR using 
appropriate primer pairs (primer sequences and PCR conditions 
described in [10]). The PCR products were sequenced by Big Dye 
Terminator cycle sequencing on an automated ABI Prism 3100 
sequencer (Applied Biosystems, Foster City, USA). 

Microarray hybridization and data normalization 

0.5-2 |J.g of total RNA from 1 1 ATC samples and from the pool 
of 19 normal thyroid tissues were labelled with Flyanin 3 or 
Hyanin 5 using the miRCURY LNA microRNA Power Lai jelling 
Kit (Exiqon), according to the manufacturer's protocol. The 
labelled RNA were purified on miRNeasy columns (Qiagen). 
Samples were hybridized using Corning Pronto! Microarray 
Hybridization Kit onto in-house-printed miRNA microarray 
slides (LNA microRNA ready-to-spot probe sets, VI 1.0, Exiqon), 
which contain 1269 miRNA (human-mouse-rat probes) spotted in 



duplicate. All hybridizations were performed with dyes swap. The 
microarray slides were washed under stringent conditions. A 
LOWESS normalization was carried out by using BRB-Array- 
Tools (Version 3.6.0). 

Search for deregulated mRNA or miRNA 

For deregulated miRNA, we have selected miRNA which 
varied from the baseline by at least 1 .5-fold in the 1 1 ATC samples 
(without any opposite regulation), called the miRNA ATC list 
(Table S4). 

For deregulated mRNA in ATC previously hybridized by our 
group, we have selected mRNA which varied from the baseline by 
at least 2-fold in all the 11 ATC samples (without any opposite 
regulation), called the mRNA ATC list [10]. 

Nonsupervised analyses 

Nonsupervised analyses were performed on the basis of 
between-sample correlation distances. Multidimensional scaling 

(MDS, as implemented by R's isoMDS function) was performed 
on all human miRNA present on the microarray. 

Immunohistochemistry 

Immunohistochemical analyses were performed on 6 ^im thick 
sections prepared from 4 paraffm-embedded ATC tissues (4 of the 
1 1 ATC studied by microarrays) and from 1 normal thyroid tissue, 
with the primary antibody directed against CD163 (Leica 
Biosystems, UK, Clone 10D6). The sections were deparaffinized, 
pretreated with CCl (EDTA, pH 8.4) and incubated with the 
antibody at a 1:100 dilution. The revelation was performed with a 
detection kit (Roche, ref 760-501). 

Validation of miRNA expression data by quantitative real- 
time RT-PCR (qRT-PCR) 

qRT-PCR was performed using TaqMan MicroRNA Assay kits 
according to the manufacturer's protocol (Apphed Biosystems). 
miRNA specific cDNA were generated from 10 ng of total RNA 
using the TaqMan microRNA RT kit and the gene-specific RT 
primers from the TaqMan microRNA Assays (Applied Biosystems) 
according to the manufacturer's protocol. Relative quantification 
was calculated based on the Pfaffl Method [18]. Ct data were 
normalized to an internal control U6 snRNA. 

AU miRNA validations were done in triplicate on 6 ATC (3 of 
the 1 1 ATC studied by microarrays and 3 independent ones). 

In situ hybridization of miRNA 

Six |im thick paraffin sections of 4 ATC tissues (4 of the 1 1 ATC 
studied by microarrays) and of 1 normal thyroid tissue were 
mounted on Super frost+glass slides and deparaffmized. The slides 
were hybridized with 40 nM of the LNA miRNA let-7g, miR-29a 
and miR-30e probes (Exiqon) and the nuclei were counterstained 
with the nuclear fast red counterstain (Vector Laboratories, 
Burlingname, CA) by Bioneer as described in [19]. A scramble 
miRNA and miR-126 were used as negative and positive controls 
respectively. LNA-miRNA sequences are described in Table S 1 . 

Search for bioinformatically predicted mRNA targets of 

miRNA 

A list of 79 mRNA targets of the 18 deregulated miRNA in 
ATC was constructed. This list resulted from the intersection of 
the bioinformatically predicted mRNA targets of the 18 deregu- 
lated miRNA (using miRTarget2 software with a miRTarget2 
score >80% [17]) and the mRNA ATC list (taking into account 
the direction of the regulation). For instance, predicted mRNA 
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Figure 1. Multidimensional scaling (MDS) of the miRNA expression data: a. from 1 1 ATC (the ATC are labeled according to their 
mutation as depicted in the figure), b. from 1 1 ATC and 5 PTC (PC), 5 AA (AA), 5 FTC (FC) and 4 FTA (FA). All human miRNA present on 
the arrays were considered for the analyses. 
doi:1 0.1 371 /journal.pone.01 03871 .gOOl 



targets of downregulated miRNA should be upregulated in our 
samples and thus were intersected with the list of mRNA 
upregulated in ATC (as described above). One of them, the 
LOX gene, predicted to be targeted by the downregulated 
miRNA-29a, was selected for the luciferase assay. 

Luciferase assay 

A PGR fragment of =400 nucleotides encompassing the miR- 
29a binding site at the end of the 3'UTR of LOX (nucleotides 
2941 to 3357) was cloned downstream of the ReniUa luciferase 
gene, in the psiCHECK2 vector (Promega) (primers used for 
plasmid construction are listed in Table S2). The HEK293 cell line 
(derived from human embryonic kidney) was used for transfection 
experiments. The cells were seeded on 24-well plates in a DMEM 
(Dulbeccos's modified eagle medium, GIBCO), L-glutamine and 
sodium pyruvate medium until 50% of confluence prior to 
transfection. 100 ng of the luciferase vector was co-transfected 
with 50 iiM of mimic miR-29a (Dharmacon) or a scramble and 
with 2 |il of LipofectAMINE 2000 (Invitrogen), according to the 
manufacturer's instructions. 24 h and 48 h after the transfection, 
luciferase activity was measured using the Dual Luciferase 
Reporter Assay system (Promega) using a MicroaLumatPlus 
(EG&G Berthold) and was normalized by measuring the firefly 
luciferase activity. Results were expressed as ratios of the 
normalized luminescence value of the cells co-transfected with 
the luciferase vector with the miR binding site wild-type or 
mutated in the seed sequence and the mimic miRNA compared to 
the normalized luminescence value of the cells co-transfected with 
the luciferase vector with the miR binding site wild-type or 



mutated in the seed sequence and the scramble. Primers used for 
directed mutagenesis are provided in Table S2. 

Statistical analyses 

TAM quantification. After Immunohistochemistry detection 
with the primary antibody directed against CD 163, 500 nucleated 
cells were counted on four different ATC tissues samples and on 
one normal adjacent tissue and the marked cells were counted for 
TAM quantification. Mean and standard deviation were then 
calculated. 

Luciferase assay. The luciferase assay was performed in 
triplicates in 3 to 4 independent cell cultures for each of the 3 
conditions (the luciferase vector with the miR binding site wild- 
type after 24 h or 48 h of transfection or the luciferase vector 
mutated in the seed sequence after 48 h of transfection). A t-test 
was performed on the ratios of normalized luminescence values (as 
described in the Luciferase assay section). 

Results 

Mutational screening of p53, BRAF, PIK3CA, H-RAS, K- 
RAS, N-RAS and p-catenin in the 1 1 ATC 

We searched for the presence of the most common mutations 
reported for ATC. On the 1 1 ATC, p53 mutation was found in 4 
(36.4%), BRAF mutation in 2 (18%), PIK3GA mutation m 1 
(10%). One sample showed both BRAF and p53 mutations 
(ATC 1). No mutation was found for RAS (H-RAS, K-RAS and N- 
RAS) or for P-catenin. 
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Figure 2. Immunohistochemistry with the CD1 63 antibody on 3 different ATC and on normal thyroid tissue, a, b, c: 20 x magnification 
showing that about 50% of the nucleated cells are TAM, d: 40 x magnification showing elongated cytoplasmic extensions in TAM as illustrated by the 
red arrows; e: 20 x magnification illustrating the almost absence of TAM in normal thyroid tissue. 
doi:1 0.1 371 /journal.pone.01 03871 .g002 



The 1 1 ATC showed distinct miRNA expression profiles 
from the other thyroid tumors 

The molecular phenotypes of ATC can best be demonstrated by 
a comprehensive microarray analysis of mRNA and miRNA 
expressions in these tissues. mRNA expression profiles of 1 1 ATC 
were obtained previously [10], and we have now performed 
miRNA expression analysis on the same tumors samples. Because 
of the absence of normal tissue counterparts for ATC, all gene 
expression profiles were compared to a common miRNA reference 
pool of 19 normal tissues which are adjacent to benign thyroid 
tumors. Overall miRNA expressions profiles from the 1 1 ATC, as 
well as from 5 PTC, 5 AA, 5 FTC and 4 FTA for comparison, 
were analyzed using a multidimensional scaling (MDS) (Figure 1). 
The MDS algorithm reduces the n-dimensions space (n: number 
of probes) into two dimensions while preserving the distances 
between the samples, and thereby visualizes the similarity 
relationships between them. MDS showed, firstly, that the ATC 
did not cluster according to their mutational status (Figure 1 a) and 
secondly, that ATC miRNA expression profiles are clearly distinct 
from those of the other thyroid tumor types (Figure lb). 

Deregulated microRNA in ATC 

ATC presented 18 commonly deregulated miRNA (Table S4), 
called the miRNA list, with 1 7 downregulated and 1 upregulated 
miRNA (selection criteria described in Material and Methods). 
These 1 8 miRNA were significantly deregulated only in the ATC 
samples, and not in the other tumor types, except for miR-144 



which was downregulated in all the thyroid tumor types (Table 
S4). 

About 60% of the nucleated cells in the ATC tissues are 
TAM 

Because TAM were shown to be present in ATC and to have a 
major role in the aggressiveness of the tumors [6], we have 
investigated their proportion by IHC with the macrophage- 
specific CD 163 antibody in 4 ATC samples studied for their 
miRNA profiles and a normal thyroid tissue as control. ATC tissue 
samples show 63%+/ — 4,6% (mean+/— standard deviation) of 
TAM among the nucleated cells (IHC are presented in 
Figure 2.a,b,c) which are almost not present in the normal thyroid 
tissue (Figure 2.e) (3%+/— 1% (mean +/— standard deviation) of 
TAM among the nucleated cells). Moreover TAM present thin 
and elongated cytoplasmic extensions forming a canopy structure 
all over the tumor tissue (as illustrated in Figure 2.d). 

Validation of miRNA expression data 

By literature. 50% of our deregulated miRNA are common 
with the miRNA modulated in the three ATC cases studied by 
Braun et al. [20]: let-7f, let7-g, miR-29a, miR-30a, miR-30e, miR- 
125b, miR-141, miR-148b and miR-200b. MiR-30a and miR- 
1 25b were also previously reported as deregulated in another study 
[21,22] (Figure SI). 

By quantitative real-time RT-PCR (qRT-PCR). To fur- 
ther validate the regulations observed by microarray, qRT-PCR 
was performed with miR-514, miR-144, miR-125b, miR-30a, let- 
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Figure 3. In situ hybridization of Let-7g, miR-29a and miR-30e in normal and ATC tissues. Tissue sections from normal and tumor thyroids 
were incubated with a full length DIG-labeled LNA probe to detect Let-7g, miR-29a and miR-30e: 20 x magnification showing a strong miRNA signal 
observed (in blue) in the normal thyrocytes and a weak signal observed in the tumor cells derived from thyrocytes. No signal was detected in the 
TAiVI. 4 ATC were investigated. Representative pictures of 3 of them are shown. 
doi:1 0.1 371 /journal.pone.01 03871 .g003 



7g, miR-29a, miR-34b, and miR-659, on 6 ATC. Three of them 
were part of the microarray tumor set, and the other three were 
new and independent ones. 

The modulations of all the miRNA were confirmed except for 
miR-34b (Figure SI). 

By in situ hybridization. Because ATC tissues are com- 
posed of about 50% of TAM which are suggested to play a major 
role in the ATC aggressiveness, we have performed in situ 
hybridization on four ATC (which were a part of the microarray 
tumor set) for 3 of the deregulated miRNA in order to localize the 
miRNA (e.g. thyrocyte derived tumor cell or TAM). 

The normal tissue section shows a strong ISH staining signal 
obtained with the LNA probes for Let-7g, miR-29a and miR-30e, 
which are located in thyrocytes. On the contrary, a very weak ISH 
signal for these 3 miRNA was observed in the four ATC samples. 
These miRNA were confirmed by a pathologist to be in tumor 
cells originated from thyrocytes and not in TAM (Figure 3). The 4 



ATC investigated showed strong staining intensities with the miR- 
126 positive control probe and almost no signal for the scramble 
probe (data not shown). 

Analysis of the putative mRNA targets of deregulated 
miRNA in the same 11 ATC samples by using DAVID 
software 

The biological action of miRNA is carried out by silencing the 
expression of their target genes. In order to highlight the role of 
the 18 miRNA which are deregulated in the ATC, a list of 79 
mRNA targets (Table S3) of the 18 deregulated miRNA (obtained 
as described in Material and Methods) was analyzed using DAVID 
software - KEGG pathway. The only significantly altered function 
(FDR score <0.05) highlighted by this analysis is called 
'cytoskeleton' and is related to the EMT process. This suggests 
that EMT, a major characteristic of ATC already observed in the 
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mRNA expression analysis [10], might be directed by the 
deregulation of miRNA in ATC, as already suggested previously 
[20]. 

Luciferase assay validates the direct interaction of LOX 
mRNA with miR-29a 

Among die mRNA targets, the LOX gene predicted to be 
targeted by the downregulated miR-29a, was selected for a 
luciferase assay in order to validate their direct interaction which 
has not been previously described. The LOX gene is a key player 
in EMT and was shown to be an interesting therapeutic target 
[23]; it is upregulated at mRNA and protein level in our ATC 
samples (data not shown). About 400 nucleotides of the 3'UTR of 
LOX containing the seed sequence of the miR-29a were cloned in 
the psiCHECK2 vector and HEK293 cells were co-transfected 
with the vector and with the mimic miR-29a or with a mimic 
scramble as control. Luciferase activity, measured after 24 h and 
48 h, was significantly decreased in the miR-29a transfected cells 
compared to the scramble transfected cells (p-value>5%). 
Deletion of a part of the seed sequence in the vector restored 
partially the luciferase activity showing the specificity of the 



binding of the miR-29a on this part of the 3'UTR of LOX 
(Figure 4). These results suggest that the miR-29a downregulates 
LOX mRNA through a direct interaction with its 3'UTR. 

Discussion 

Many studies support the idea that miRNA play a role in 
tumorigenesis. This study was undertaken to further characterize 
the molecular mechanisms underlying ATC pathogenesis. To 
achieve this goal, the same 1 1 ATC samples previously analyzed 
for mRNA expression [10] have now been investigated for 
miRNA expression allowing a direct correlation between both 
expression profiles, as each miRNA can regulate many mRNA. 
MiRNA are often considered as fine-tunes regulators of the 
genome. 

miRNA expression analysis was performed on 1 1 ATC and on 
19 other less aggressive thyroid tumors for comparison including 
PTC, FTC, FA and AA. 1 7 downregulated miRNA and only one 
upregulated miRNA were identified as commonly deregulated in 
all ATC compared to a normal thyroid pool. This general 
downregulation of miRNA is also observed in many other cancers 
[24]. While the patterns of mRNA expressions in PTC and ATC 
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Figure 5. Proposed scheme synthetizing the miRNA and mRNA expression analyses of the same 1 1 ATC and highlighting the EMT 
feature of this tumor. Only direct binding of miRNA on mRNA are represented (green: downregulated m(i)RNA; red: upregulated m(i)RNA. 
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are largely overlapping and show a distinct, but mostly quantita- 
tive evolution from one to the other [10], the overall downreg- 
ulation of miRNA in ATC and their specific miRNA expression 
not found in the other tumor types can better explain the dramatic 
qualitative switches in pathology and clinical evolution. Another 
hypothesis explaining this physiopathology switch might be the 
high presence of TAM observed by others [6] and confirmed by 
immunohistochemistry in our ATC tissues (about 60% of the 
nucleated cells), while the other tumor types are less infiltrated by 
TAM [6] . TAM show thin, elongated cytoplasmic extensions that 
form a canopy structure aU over timior tissue. Paracrine or direct 
communication via the cytoplasmic extensions between the TAM 
and the thyroid tumor cells are often mentioned and are suggested 
to enhance tumor progression [25-27]. In situ hybridization 
experiments have assessed that the downregulated let-7g, miR-29a 
and miR-30e are present in tumor cells derived from thyrocytes 
and not in TAM. In the case of very heterogeneous tumors, such 
as ATC, it is thus essential to carefully verify the exact cellular 
localization of any transcript or protein identified by whole tissues 
analyses, before drawing any conclusion on its physiological role. 

Studying the function of the mRNA targets of the deregulated 
miRNA should highlight a role of these miRNA in the 
pathogenesis of the ATC. The KEGG analysis pointed out the 
EMT process. Indeed, ATC are highly invasive, and part of 
patients die because of overwhelming distant metastases disease 
involving an EMT process. Moreover, the ATC tissue is 
completely dedifferentiated, without any structure anymore, and 
the cells present a mesenchymal morphology. The EMT process is 



the consequence of a molecular reprogramming of the cells 
consisting in an increased expression of mesenchymal type genes, 
like the mesenchymal N-cadherin (CDH2) which facilitates tumor 
cell binding to the stroma, and in a decreased expression of 
epithelial type genes, like E-cadherin (CDHl) which allows cell- 
cell interactions. TAM can induce the EMT process for example 
through paracrine TGF-(3 signalling and consecutive activation of 
the fi-catenin pathway [26] or through TLR4/IL10 signaling [25]. 
The EMT switch may also be the consequence of the action of 
various miRNA modulating the expression of different mRNA, as 
synthetized in Figure 5. 

In the following paragraphs, miRNA regulations which are not 
significant in all the ATC are specified by the fraction of positive 
tumors (no specification for those significantly alterated in all our 
1 1 ATC). In our mRNA data [10], CDHl is downregxilated and 
CDH2 is upegulated in ATC [10]. CDHl is transcriptionally 
repressed by SNAIL, ZEBl, ZEB2 and TWIST, and induces 
subsequent breakdown of cell junction and the ensuing loss of ceU 
polarity [28]. SNAIL, ZEBl, ZEB2 and TWIST are all 
upregulated in ATC samples [10]. Three members of the miR- 
200 family, miR-141, miR-200b and miR-429, downregulated in 
our ATC (11/11, 11/11 and 9/11 respectively), have been shown 
to silence the ZEB proteins [29-32]. In addition, miR-101 
(downregulated in 10/11 ATC) targets EZH2 (upregulated in all 
the ATC [33]), and mediates transcriptional silencing of CDHl 
[34]. 

Activation of the RAS/MAPK pathway also induces EMT, 
together with HMGA2 upregulation. Let-7, for which many 
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members of this family of miR are downregulated in our ATC 
samples, inversely regulates HMGA2 [35], upregulated in our data 
at mRNA and protein levels (data not shown). HMGA2 itself 
induces the expression of SNAIL, a transcriptional repressor of E- 
cadherin [35,36]. This procxss can be initiated, for example, via 
the oncostatin M which can be secreted by M2 type macrophages 
and which downregulates miR-200 and let-7 family members and 
finally induces the EMT process [36] . Also, mir-98, downregulated 
in all our ATC tumors, negatively regulates HMGA2. The 
reintroduction of miR-98 accelerated the inhibition of glioma cell 
invasion and can therefore be interesting for furth("r characteriza- 
tion for ATC therapy [37]. Another miRNA promoting the EMT 
process is miR-335 (downregulated in 10/11 ATC) which targets 
SOX4 and TNC (upregulated in 10/11 and 11/11 ATC 
respectively) [38,39]. Moreover, we have now validated the direct 
interaction between the downregulated miR-29a and the upregu- 
lated LOX mRNA, which is also an EMT key player [23], 
suggesting another way for cancer cells to promote EMT. 

Benefits obtained from chemotherapy and radiation therapy of 
ATC remain only marginal and there are no alternative 
treatments yet [4,5], therefore new therapeutic approaches are 
needed. Our results suggest several approaches: first, detection of a 
circulating miRNA in the plasma, present in a remarkably stable 
form because protected from endogenous RNase activity, may 
allow to predict diagnosis or prognosis and eventually to treat the 
patients earlier [40] . This could be tested for the upregulated miR- 
659. Second, reinduction of one or more miRNA, downregulated 
in tumor cells, may be an interesting approach to treat the ATC 
patients but, of course, it raises the question of cell specificity. 
Third, as our results highlight a high proportion of TATvI in ATC, 
and as TAM is pro-tumorigenic in advanced thyroid cancer, they 
could be then targeted pharmacologically [41]. 

Supporting Information 

Figure SI Expression levels of 13 miRNA measured by qRT- 
PCR in 3 ATC studied by microarrays and in 3 ATC independant 

References 

1. Smallridge RG, Marlow LA, Copland JA (2009) Anaplastic thyroid cancer: 
molecular padiogenesis and emerging dicrapics. Endocr Relat Cancer 16: 17— 
44, 

2. Ain KB (1999) Anaplastic thyroid carcinoma: a therapeutic challenge. Semin 
Surg Oncol 16: 64-69. 

3. Nagaiah G, Hossain A, Mooncy CJ, Parmcnticr J, Rcmick SC (201 1) Anaplastic 
thyroid cancer: a review of epidemiology, pathogenesis, and treatment. J Oncol 
2011: 542358. 

4. Are C, Shaha .\R (2006) .\naplastie thyroid carcinoma: biology, pathogenesis, 
prognostic laetors, and treatment approaches. .Ann Surg Oneol 13: 453—464. 

5. Cooper DS, Doherty GM, Haugen BR, Kloos RT, Lee SL, et al. (2006) 
Management guidelines for patients with thyroid nodules and differentiated 
thyroid cancer. Thyroid 16: 109-142. 

6. Caillou B, Talbot M, Weyemi U, Koche-Durieu C, Al OA, et al. (201 1) Tumor- 
associated macrophages (TAMs) form an interconnected cellular supportive 
network in anaplastic thyroid carcinoma. PLoS One 6: e22567. 10.1371/ 
journal.pone.0022567 [doi];PONE-D- 10-05881 [pii]. 

7. Kim S, Gho S\V, Min HS, Kim KM, Yeom GJ, et al. (2013) The expression of 
tumor-assoeiated macrophages in papillary thyroid carcinoma. Endocrinol 
Metab (Seoul) 28: 192-198. 'l0.3803/EnM.2013. 28. 3. 192 [doij. 

8. Ryder M, (ihossein RA, Riearte-Filho JC, KnaufJA, EaginJ.A (2008) Increased 
density of tumor-associated macrophages is associated with decreased survival in 
advanced thyroid cancer. Endocr Relat Cancer 15: 1069-1074. ERC-08-0036 
|pii];10.1677/ERC-08-0036 [doi]. 

9. Qing W, Fang WY, Ye L, Shen LY, Zhang XF, et al. (2012) Density of tumor- 
associated macrophages correlates with lymph node metastasis in papillary 
tiiyroid carcinoma. Thyroid 22: 905-910. 10.1089/thy.201 1.0452 [doi]. 

10. Hebrant A, Dom G, Dewaele M, Andry G, Tresallet C, et al. (2012) mRNA 
expression in papillary and anaplastic thyroid carcinoma: molecular anatomy of 
a killing switch. PLoS One 7: e37807. 

11. He L, Hannon GJ (2004) MicroRNAs: small RNAs with a big role in gene 
regulation. Nat Rev Genet 5: 522-531. 10.1038/nrgl379 [doi];nrgl379 [pii]. 



or obtained fi-om the literature. The microarray expressions are 
included for comparision. Log2 ratios represent the expression 
ratios of the genes in the tumors versus the normal tissues. 
(PPT) 

Table SI 

(XLSX) 

Table S2 

PCLS) 

Table S3 

(XLSX) 

Table S4 Deregulated miRNA in tbe 11 ATC samples 
and their corresponding values in papillary carcinoma 
(PC), follicular carcinoma (FC), follicular adenoma (FA) 
and autonomous adenoma (AA). MiRNA expression values 
are in bold if they varied from the baseline by at least 1.5-fold (log2 
of expression ratios). 
(PPTX) 

Acicnowledgments 

We thank Chaiilal Dcgracl" lor excellent technical assistance and Jacques 
Emile DumonL lor helplul and eonstrtictive discussions. 

ATC tissue samples were obtained from the tissue bank of Regional 
Reference Cancer Center of Lille "Tumorotheqne du centre regional de 
Reference en Canct^rologie" (France), of the Bordet Institute (Brussels, 
Belgium) and of the Cancer Center of the Cliniques Universitaires Saint- 
Luc (UCL, Brussels, Belgium), project CDCUGLR. 

Author Contributions 

Conceived and designed the experiments: AH CM. Performed the 
experiments: AH AD BS. Analyzed the data: AH. Contributed reagents/ 
materials/analysis tools: SF MS AA CL NdSA EL GA. Wrote the paper: 
AH. 



12. Kiriakidou M, Nelson PT, Kouranov A, Fitziev P, Bouyioukos C, et al. (2004) A 
combined computational-experimental approach predicts human microRNA 
targets. Genes Dev 18; 1165-1178. 

13. Ba^a S, Bracht J, Hunter S, Massirer K, HoltzJ, et al. (2005) Regulation by let- 
7 and lin-4 miRNAs results in target mRNA degradation. Cell 122: 553-563. 

14. Du T, Zamore PD (2005) microPrimcr: the biogenesis and fimction of 
microRNA. Development 132: 4645-4652. 

15. Djuranovic S, Nahvi A, Green R (2011) A parsimonious model for gene 
regulation by miRNAs. Science 331: 550-553. 

16. Guo H, Ingolia NT, VVeissman JS, Bartel DP (2010) Mammalian mieroRNAs 
predominantly act to decrease target mRNA levels. Nature 466: 835-840. 

17. Wang X, El Naqa IM (2008) Prediction of both eonser\ed and nonconserved 
microRNA targets in animals. Bioinformatics 24: 325-332. 

18. PfafQ MW (2001) A new mathematical model for relative quantification in real- 
time RT-PCR. Nucleic Acids Res 29: e45. 

19. Nielsen BS, Jorgensen S, FogJU, Sokilde R, Christensen IJ, et al. (2011) High 
levels of microRNA-2 1 in the stroma of colorectal cancers predict short disease- 
free sur\?ival in stage 11 colon cancer patients. Clin Exp Metastasis 28: 27—38. 
1 0. 1 007 /s 1 0585-0 1 0-9355-7 [doi] . 

20. Braun J, Hoang-Vu C, Dralle H, Huttelmaier S (2010) Downregulation of 
microRN.As directs the EMT and invasive potential of anaplastic thyroid 
carcinomas. Oncogene 29: 4237-4244. 

21. Visone R, Pallante P, Vecchione .A, Cirombella R, Ferraein M, et al. (2007) 
Specific microRNAs are downregulated in human thyroid anaplastic carcino- 
mas. Oncogene 26: 7590-7595. 

22. Pallante P, Battista S, Pierantoni GM, Fusco A (2014) Deregulation of 
microRNA expression in thyroid neoplasias. Nat Rev Endocrinol 10: 88-101. 
nrendo.2013.223 [pii];10.1038/nrendo.2013.223 [doi]. 

23. Sehietke R, Warnecke C, Wacker 1, SchodelJ, Mole DR, et al. (2010) The lysyl 
oxidases LOX and LOXL2 are necessary and sufficient to repress E-eadherin in 
hypoxia: insights into cellular transformation processes mediated by HIF-1. 



PLOS ONE I www.plosone.org 



8 



August 2014 I Volume 9 | Issue 8 | e103871 



miRNA Expression in ATC 



J Biol Chem 285: 6658-6669. M109.042424 [^ii];10.1074/jbc.M109.042424 
[doi]. 

24. Hammond SM (2006) MicroRNAs as oncogenes. Curr Opin Genet Dev 16: 4- 
9. 

25. liu CY, XuJY, Shi XY, Huang W, Ruan TY, et al. (2013) M2-polarized tumor- 
associated macrophages promoted epithelial-mesenchymal transition in pancre- 
atic cancer cells, partially through TLR4/IL-10 signaling pathwav. Lab Invest 
93: 844-854. Ialiuur^l2013f)9 [pii]:10. 1038/labimcsl. 201 3.09 [doij. 

26. Bondc AK, Tischlcr \', Kumar S, Soltcrmann A, Schwcndcncr Rj\ (2012) 
Intratumoral macrophages contribute to epithelial-mesenchymal transition in 
solid tumors. BMC Cancer 12: 35. 1471-2407-12-35 [piiJ;lo'.1186/1471-2407- 
12-35 [doi]. 

27. Baay M, Brouwer A, Pauwels P, Peeters M, Lardon F (2011) Tumor cells and 

tumor-associated macrophages: secreted proteins as potential targets for therapy. 
Clin Dev Immunol 2011: ,")65187. 10.1 155/20Il/,lf),1187 [doi/. 

28. Yang J, Weinberg RA (2008) Epithelial-mesenchymal transition: at the 
crossroads of development and tumor metastasis. Dev Cell 14: 818-829. 

29. Aigner A (2011) MicroRNAs (miRNAs) in cancer invasion and metastasis: 
therapeutic approaches based on metastasis-related miRXAs. J Alol Aled. 

30. Dykxhoorn DM, VVu Y, Xie H, Yu F, Lai A, et al. (2009) miR-200 enhances 
mouse breast cancer cell colonization to form distant metastases. PLoS One 4: 
e7181. 

31. Gregory PA, Bert AG, Paterson EL, Barry SC, Tsykin A, et al. (2008) The miR- 
200 family and miR-205 regulate epithelial to mesenchymal transition by 
targeting ZEBl and SIPl. Nat CeU Biol 10: 593-601. 

32. Saycd D, Abdcllatif M (2011) MicroRNAs in development and disease. Physiol 

Rev 91: 827-887. 

33. Varambally S, Cao Q, Mani RS, Shankar S, Wang X, et al. (2008) Genomic loss 
of microRNA-101 leads to overexpression of histone methyltransferase EZH2 in 
cancer. Science 322: 1695-1699. 



34. Cao Q, Yu J, Dhanasekaran SM, Kim JH, Mam RS, et al.(2008) Repression of 
E-cadherin by the polycomb group protein EZH2 in cancer. Oncogene 27: 
7274-7284. 

35. Watanabe S, Ueda Y, Akaboshi S, Hino Y, Sekita Y, et al. (2009) HMGA2 
maintains oncogenic RAS-induced epithelial-mesenchymal transition in human 
pancreatic cancer cells. Am J Pathol 174: 854—868. 

36. Guo L, Chen C, Shi M, Wang F, Chen X, et al.(201 3) Stat3-coordinatcd Lin-28- 
let-7-HMGA2 and miR-200-ZEBl circuits initiate and maintain oncostatin M- 
driven cpithelial-mcscnchvmal transition. Oncogene 32: 5272—5282. 
onc2012573 [piij;10.1038/onc.2012.573 [doiJ. 

37. Chen Z, Cheng Q, Ma Z, Xi H, Peng R, et al. (2013) Overexpression of RKIP 
Inhibits Cell Invasion in Glioma Cell Lines through Uprcgulation of miR-98. 
Biomed Res Int 2013: 695179. 10. 1 155/201 3/69,')l 79 [doi]. 

38. Tavazoie SF, Alarcon C, Oskarsson T, Paduii D, Wang Q, et al. (2008) 
Endogenous human microRNAs that suppress breast cancer metastasis. Nature 
451: 147-152. namre06487 [pii];10.1038/nature06487 [doi]. 

39. Pug KJ, Yoshida M, Zhang XH, Shu W, Lee H, et al. (201 1) MieroRNA-335 
inhibits tumor reinitiation and is silenced through genetic and epigenetic 
mechanisms in human breast cancer. Genes Dev 25: 226-231. 25/3/226 
[pii];I0.II01/gad.I974211 [doi]. 

40. MitcheU PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, et al. (2008) 
Circulating microRNAs as stable blood-based markers for cancer detection. Proc 
Natl Acad Sci USA 105: 1051,3-10518. 

41. Ryder M, Gild M, Hohl TM, Pamer E, Knauf J, et al. (2013) Genetic and 
pharmacological targeting of CSF-l/CSF-lR inhibits tumor-associated macro- 
phages and impairs BRAF-induced thyroid cancer progression. PLoS One 8: 
e54302. 10.I37I/joumal.pone,0054302 [doi] ;PONE-D- 12-29549 [pii]. 



PLOS ONE I www.plosone.org 



9 



August 2014 I Volume 9 | Issue 8 | e103871 



